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Abstract. We study the time evolution of ranking and spectral properties of the Google matrix of English 
Wikipedia hyperlink network during years 2003 - 2011. The statistical properties of ranking of Wikipedia 
articles via PageRank and CheiRank probabilities, as well as the matrix spectrum, are shown to be stabi- 
lized for 2007 - 2011. A special emphasis is done on ranking of Wikipedia personalities and universities. We 
show that PageRank selection is dominated by politicians while 2DRank, which combines PageRank and 
CheiRank, gives more accent on personalities of arts. The Wikipedia PageRank of universities recovers 80 
percents of top universities of Shanghai ranking during the considered time period. 
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PACS. 89.75.Fb Structures and organization in complex systems 
trees - 89.20.Hh World Wide Web, Internet 



89. 75. He Networks and genealogical 



1 Introduction 

At present Wikipedia [T] became the world largest En- 
cyclopedia with open public access to its contain. A re- 
cent review [2J represents a detailed description of pub- 
lications and scientific research of this modern Library 
of Babel, which stores an enormous amount of informa- 
tion, approaching the one described by Jorge Luis Borges 
[3] . The hyperlinks of citations between Wikipedia articles 
represent a directed network which reminds the structure 
of the World Wide Web (WWW). Hence, the mathemat- 
ical tools developed for WWW search engines, based on 
the Markov chains [4| , Perron- Frobenius operators 5^ and 
the PageRank algorithm of the corresponding Google ma- 
trix [6l[7|, give solid mathematical grounds for analysis of 
information flow on the Wikipedia network. In this work 
we perform the Google matrix analysis of Wikipedia net- 
work of English articles extending the results presented 
in [51l^, |10l[TT| . The main new element of this work is 
the study of time evolution of Wikipedia network during 
the years 2003 to 2011. We analyze how the ranking of 
Wikipedia articles and the spectrum of the Google matrix 
G of Wikipedia are changed during this period. 

The directed network of Wikipedia articles is const- 
ructed in a usual way: a directed link is formed from an 
article j to an article i when j quotes i and an element Aij 
of the adjacency matrix is taken to be unity when there is 
such a link and zero in absence of link. Then the matrix 
Sij of Markov transitions is constructed by normalizing 
elements of each column to unity {J2j ^ij ^1) ^^^ re- 
placing columns with only zero elements {dangling nodes) 



by 1/iV, with N being the matrix size. Then the Google 
matrix of the network takes the form OIT] : 



G,y = aS^j + (1 - a)/N 



(1) 



The damping parameter a in the WWW context describes 
the probability (1 — a) to jump to any node for a random 
surfer. For WWW the Google search engine uses a ~ 0.85 
[7j . The matrix G belongs to the class of Perron- Frobenius 
operators [SJ[7], its largest eigenvalue is A = 1 and other 
eigenvalues have |A| < a. The right eigenvector at A = 1, 
which is called the PageRank, has real nonnegative ele- 
ments P{i) and gives a probability P{i) to find a random 
surfer at site i. It is possible to rank all nodes in a de- 
creasing order of PageRank probability P{K{i)) so that 
the PageRank index K{i) counts all N nodes i according 
their ranking, placing the most popular articles or nodes 
at the top values K = 1,2, 3.... 

Due to the gap 1 — a « 0.15 between the largest eigen- 
value A = 1 and other eigenvalues the PageRank algo- 
rithm permits an efficient and simple determination of the 
PageRank by the power iteration method [7]. It is also 
possible to use the powerful Arnoldi method [l2l[T3] . [I4] 
to compute efficiently the eigenspectrum Ai of the Google 
matrix: 



N 



^Gjki^iik) = Xiiptij) 



(2) 



fc=i 



The Arnoldi method allows to find a several thousands 
of eigenvalues Xi with maximal |A| for a matrix size N 
as large as a few tens of millions [TU1ITT| , [T^[T5] . Usually, 
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at a = 1 the largest eigenvalue A = 1 is highly degen- 
erate |15j due to many invariant subspaces which define 
many independent Perron-Frobenius operators providing 
(at least) one eigenvalue A = 1. 

In addition to a given directed network Aij it is use- 
ful to analyze an inverse network with inverted direction 
of links with elements of adjacency matrix Aij — > Aji. 
The Google matrix G* of the inverse network is then con- 
structed via corresponding matrix S* according to the re- 
lations ([T]) using the same value of a as for the G matrix. 
This time inversion approach was used in [16,17 but the 
statistical properties and correlations between direct and 
inversed ranking were not analyzed there. In [TS], on an 
example of the Linux Kernel network, it was shown thus 
this approach allows to obtain an additional interesting 
characterization of information fiow on directed networks. 
Indeed, the right eigenvector of G* at eigenvalue A = 1 
gives a probability P*{i), called CheiRank vector [8 . It 
determines a complementary rank index K*{i) of network 
nodes in a decreasing order of probability P*{K*{i)) [51 
[5]. |10l[T5] . It is known that the PageRank probability is 
proportional to the number of ingoing links characteriz- 
ing how popular or known is a given node. In a similar 
way the CheiRank probability is proportional to the num- 
ber of outgoing links highlighting the node communicativ- 
ity (see e.g. [71ITO]. pH l l ^ . piS] ) . The statistical properties 
of distribution of indexes K{i),K*{i) on the PageRank- 
CheiRank plane are described in [9]. 

In this work we apply the above mathematical meth- 
ods to the analysis of time evolution of Wikipedia network 
ranking using English Wikipedia snapshots dated by De- 
cember 31 of years 2003, 2005, 2007, 2009, 2011. In ad- 
dition we use the snapshot of August 2009 (200908) ana- 
lyzed in [S] . The parameters of networks with the number 
of articles (nodes) N, number of links Ni and other in- 
formation are given in Tables 1,2 with the description of 
notations given in Appendix. 

The paper is composed as following: the statistical 
properties of PageRank and CheiRank are analyzed in 
Section 2, ranking of Wikipedia personalities and univer- 
sities are considered in Sections 3, 4 respectively, the prop- 
erties of spectrum of Google matrix are considered in Sec- 
tion 5, the discussion of the results is presented in Section 
6, Appendix Section 7 gives network parameters. 



2 CheiRank versus PageRank 

The dependencies of PageRank and CheiRank probabili- 
ties P{K) and P*{K*) on their indexes K, K* at different 
years are shown in Fig. [T] The top positions of K are oc- 
cupied by countries starting from United States while at 
the top positions of K* we find various listings (e.g. ge- 
ographical names, prime ministers etc.; in 2011 we have 
appearance of listings of listings). Indeed, the countries 
accumulate links from all types of human activities and 
nature, that make them most popular Wikipedia articles, 
while listings have the largest number of outgoing links 
making them the most communicative articles. 




10° 10' 10^^ 10^ 10* 10^ lO'^ 10' 10° 10^ 10^ 10^ 10* 10^ 10*^ io' 



Fig. 1. PageRank probability P{K) (left panel) and CheiRank 
probability P''{K*) (right panel) are shown as a function of the 
corresponding rank indexes K and K* for English Wikipedia 
articles at years 2003, 2005, 2007, 200908, 2009, 2011; here the 
damping factor is a = 0.85. 



The data of Fig. [T] show that the global behavior of 
P{K) remains stable from 2007 to 2011. The probability 
P*{K*) is stable in the time interval 2007 - 2009 while at 
2011 we see the appearance of peak at 1 < K* < 10 that 
is related to introduction of listings of listings which were 
absent at earlier years. At the same time the behavior of 
P*{K*) in the range IQ < K* < 10^ remains stable for 
2007- 2011. 

Each article i has its PageRank and CheiRank indexes 
K(i)^ K*{i) so that all articles are distributed on two- 
dimensional plane of PageRank-CheiRank indexes. Fol- 
lowing [Slin] we present the density of articles in the 2D 
plane {K, K*) in Fig. [2] The density is computed for 100 x 
100 logarithmically equidistant cells which cover the whole 
plane {K,K*) for each year. The density distribution is 
globally stable for years 2007-2011 even if there are arti- 
cles which change their location in 2D plane. We see an 
appearance of a mountain like ridge of probability along a 
line \b.K* Ki \nK + 4.6 that indicate the presence of cor- 
relation between P{K{i)) and P*{K*{i)). Following ^^ 
[T5] we characterize the interdependence of PageRank and 
CheiRank vectors by the correlator 



N 



NY,P{K{^))P*{K*{^)) 



(3) 



We find the following values of the correlator at vari- 
ous time slots: n = 2.837(2003), 3.894(2005), 4.121(2007), 
4.084(200908), 6.629(2009), 5.391(2011). During that pe- 
riod the size of the network increased almost by 10 times 
while K, increased less than 2 times. This confirms the sta- 
bility of the correlator k, during the time evolution of the 
Wikipedia network. 

In the next two Sections we analyze the time variation 
of ranking of personalities and universities. 



3 Ranking of personalities 

To analyze the time evolution of ranking of Wikipedia 
personalities (persons or humans) we chose the top 100 
persons appearing in the ranking list of Wikipedia 200908 
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Fig. 2. Density of Wikipedia articles in the CheiRank ver- 
sus PageRank plane at different years. Color is proportional to 
logarithm of density changing from minimal nonzero density 
(dark) to maximal one (white), zero density is shown by black 
(distribution is computed for 100 x 100 cells equidistant in log- 
arithmic scale; bar shows color variation of natural logarithm 
of density); left column panels are for years 2003, 2007, 200908 
and right column panels are for 2005, 2009, 2011 (from top to 
bottom). 
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given in [5J in order of PageRank, CheiRank and 2DRank. 
We remind that 2DRank K2 is obtained by counting nodes 
in order of their appearance on ribs of squares in (K, K*) 
plane with their size growing from K — 1 to K = N [S]. 

The distributions of personahties in PageRank-CheiRank 
plane is shown at various time slots in Fig. [3] There are 
visible fluctuations of distribution of nodes for years 2003, 
2005 when the Wikipedia size has rapid growth. For other 
years the distribution of top 100 nodes of PageRank and 
2DRank is stable even if individual nodes change their 
ranking. For top 100 of CheiRank the fluctuations remain 
strong during all years. Indeed, the number of outgoing 
links is more easy to be modified by authors writing a 
given article, while a modification of ingoing links depends 
on authors of other articles. 

In Fig. [3] we also show the distribution of top 100 per- 
sonalities from Hart's book [32] (the list of names is also 
available at the web page [8]). This distribution also re- 
mains stable in years 2007-2011. It is interesting to note 
that while top PageRank and 2DRank nodes form a kind 



Fig. 3. Change of locations of top-rank persons of Wikipedia 
in K-K* plane. Each list of top ranks is determined by data 
of top 100 personalities of time slot 200908 in corresponding 
rank. Data sets are shown for (a) PageRank, (b) CheiRank, (c) 
2DRank, (d) rank from Hart g2]. 



of droplet in {K, K*) plane, the distribution of Hart's per- 
sonalities approximately follows the ridge along the line 
InK* ^lnK + A.6. 

The time evolution of top 10 personalities of slot 200908 
is shown in Fig. U for PageRank and 2DRank. For PageR- 
ank the main part of personalities keeps their rank posi- 
tion in time, e.g. G.W.Bush remains at first-second po- 
sition. B.Obama significantly improves his ranking as a 
result of president elections. There are strong variations 
for Elizabeth II which we relate to modification of arti- 
cle name during the considered time interval. We also see 
a steady improvement of ranking of C.Linnaeus that we 
attribute to a growth of various botanic descriptions and 
listings at Wikipedia articles which quote his name. For 
2DRank we observe stronger variations of K2 index with 
time. Such a politician as R.Nixon has increasing K2 in- 
dex with time since the period of his presidency goes in 
the past. At the same time singers and artists remain at 
approximately constant level of K2 ■ 

In |5] it was pointed out that the top personalities of 
PageRank are dominated by politicians while for 2DRank 
the dominant component of human activity is represented 
by artists. We analyze the time evolution of the distri- 
bution of top 30 personalities over 6 categories of human 
activity (politics, arts, science, religion, sport and etc (or 
others)). The category etc contains only C.Columbus. The 
results are presented in Fig. [SJ They clearly show that 
the PageRank personalities are dominated by politicians 
whose percentage increases with time, while the percent 
of arts decreases. For 2DRank we see that the arts are 
dominant even if their percentage decreases with time. 
We also see the appearance of sport which is absent in 
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Fig. 4. Time evolution of top 10 personalities of year 200908 
in indexes of PageRank K (a) and 2DRank A'2 (b); B.Obama 
is added in panel (a). 



PageRank. The mechanism of the qualitative ranking dif- 
ferences between two ranks is related to the fact that 
2DRank takes into account via CheiRank a contribution 
of outgoing links. Due to that singers, actors, sportsmen 
increase their ranking since they are listed in various mu- 
sic albums, movies sport competition results. Due to that 
the component of arts gets higher positions in 2DRank in 
contrast to politics dominance in PageRank. Thus the two- 
dimensional ranking on PageRank-CheiRank plane allows 
to select qualities of nodes according to their popularity 
and communicativity. 



4 Ranking of universities 

The local ranking of top 100 universities is shown in Fig. IH] 
for years 2003, 2005, 2007 and in Fig.[7]for 2009, 200908, 
2011. The local ranking is obtained by selecting top 100 
universities appearing in PageRank listing so that they get 
their university ranking K from 1 to 100. The same proce- 
dure is done for CheiRank listing of universities obtaining 
their local CheiRank index K* from 1 to 100. Those uni- 
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Fig. 5. Left panel: distribution of top 30 PageRank personal- 
ities over 6 activity categories at various years of Wikipedia. 
Right panel: distribution of top 30 2DRank personalities over 
the same activity categories at same years. Categories are pol- 
itics, art, science, religion, sport, etc (other). Color shows the 
number of personalities for each activity expressed in percents. 



versities which enter inside 100 x 100 square on the local 
index plane {K, K*) are shown in Figs. [HI El 

The data show that the top PageRank universities are 
rather stable in time, e.g. U Harvard is always on the first 
top position. At the same time the positions in K* are 
strongly changing in time. To understand the origin of this 
variations in CheiRank we consider the case of U Cam- 
bridge. Its Wikipedia article in 2003 is rather short but it 
contains the list of all 31 Colleges with direct links to their 
corresponding articles. This leads to a high position of U 
Cambridge with university K* — 4 in 2003 (Fig. |8]). How- 
ever, with time the direct links remain only to about 10 
Colleges while the whole number of Colleges are presented 
by a list of names without links. This leads to a significant 
increase of index up to K* k, 40 at Dec 2009. However, 
at Dec 2011 U Cambridge again improves significantly its 
CheiRank obtaining K* — 2. The main reason of that 
is the appearance of section of "Notable alumni and aca- 
demics" which provides direct links to articles about out- 
standing scientists studied and/or worked at U Cambridge 
that leads to second position at K* = 2 among all uni- 
versities. We note that in 2011 the top CheiRank Univer- 
sity is George Mason University with university K* = 1. 
The main reason of this high ranking is the presence of 
detailed listings of alumni in politics, media, sport with 
direct links to articles about corresponding personalities 
(including former director of CIA). These two examples 
show that the links, kept with a large number of univer- 
sity alumni, significantly increase CheiRank position of 
university. We note that artistic and politically oriented 
universities usually preserve more links with their alumni. 

The time evolution of global ranking of top 10 univer- 
sities of year 200908 for PageRank and 2DRank is shown 
in Fig. |S1 The results show the stability of PageRank or- 
der with a clear tendency of top universities (e.g. Harvard) 
to go with time to higher and higher top positions of K . 
Thus for U Harvard the global value of K changes from 
K « 300 in 2003 to K k. 100 in 2011, while the whole size 
N of the Wikipedia network increases almost by a fac- 
tor 10 during this time interval. Since Wikipedia ranks all 
human knowledge, the stable improvement of PageRank 
indexes of universities reflects the global growing inipor- 
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Fig. 6. University of Wikipedia articles in the local CheiRank 
versus PageRank plane at different years; panels are for years 
2003, 2005, 2007 (from top to bottom). 
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5 Google matrix spectrum 

Finally we discuss the time evolution of the spectrum of 
Wikipedia Google matrix taken at a = 1. We perform the 
numerical diagonalization based on the Arnoldi method 
[121113) using the additional improvements described in 
[HlfTS] with the Arnold dimension ua = 6000. The Google 
matrix is reduced to the form 
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Fig. 8. Time evolution of global ranking of top 10 Universities 
of year 200908 in indexes of PageRank K (a) and 2DRank K2 
(b). 



tance of universities in the world of human activity and 
knowledge. 

The time evolution of the same universities in 2DRank 
remains stable in time showing certain interchange of their 
ranking order. We think that an example of U Cambridge 
considered above explains the main reasons of these fluc- 
tuations. In view of 10 times increase of the whole network 
size during the period 2003 - 2011 the average stability of 
2DRank of universities also confirms the significant im- 
portance of their place in human activity. 

Finally we compare the Wikipedia ranking of universi- 
ties in their local PageRank index K with those of Shang- 
hai university ranking |23] . In the top 10 of Shanghai uni- 
versity rank the Wikipedia PageRank recovers 9 (2003), 9 
(2005), 8 (2007), 7 (2009), 7 (2011). This shows that the 
Wikipedia ranking of universities gives the results being 
very close to the real situation. A small decrease of overlap 
with time can be attributed to earlier launched activity of 
leading universities on Wikipedia. 



5*, 
S 



ss '-^sc 



(4) 



where Sgs describes disjoint subspaces Vj of dimension dj 
invariant by applications of S; Sec depicts the remaining 
part of nodes forming the wholly connected core space. 
We note that Sss is by itself composed of many small di- 
agonal blocks for each invariant subspace and hence those 
eigenvalues can be efficiently obtained by direct ( "exact" ) 
numerical diagonalization. The total subspace size iV^, the 
number of independent subspaces Nd, the maximal sub- 
space dimension dmax and the number Ni of S eigenvalues 
with A = 1 are given in Table 2 (See also Appendix). The 
spectrum and eigenstates of the core space Sec are de- 
termined by the Arnoldi method with Arnoldi dimension 
UA giving the eigenvalues Xi of Sec with largest modulus. 
Here we restrict ourselves to the statistical analysis of the 
spectrum A^. The analysis of eigenstates ij^i {Gipi = Xiipi), 
which has been done in |TT] for the slot 200908, is left for 
future studies. 

The spectrum for all Wikipedia time slots is shown 
in Fig. M for G and in Fig. [TO] for G* . We see that the 
spectrum remains stable for the period 2007 - 2001 even if 
there is a small difference of slot 200908 due to a slightly 
different cleaning link procedure (see Appendix). For the 
spectrum of G* in 2007 - 2001 we observe a well pro- 
nounced 3-6 arrow star structure. This structure is very 
similar to those found in random unistochastic matrices of 
side 3-4 [53] (see Fig. 4 therein) . This fact has been pointed 
in [TTl for the slot 200908. Now we see that this is a generic 
phenomenon which remains stable in time. This indicates 
that there are dominant groups of 3-4 nodes which have 
structure similar to random unistochastic matrices with 
strong ties between 3-4 nodes and various random permu- 
tations with random hidden complex phases. The spectral 
arrow star structure is significantly more pronounce for 
the case of G* matrix. We attribute this to more signif- 
icant fluctuations of outgoing links that probably makes 
sectors of G* to be more similar to elements of unistochas- 
tic matrices. A further detailed analysis will be useful to 
understand these arrow star structure and its links with 
various communities inside Wikipedia. 

As it is shown in [TT the eigenstates of G and G* 
select certain well defined communities of the Wikipedia 
network. Such an eigenvector detection of the communi- 
ties provides a new method of communities detection in 
addition to more standard methods developed in network 
science and described in [3S]. However, the analysis of 
eigenvectors represents a separate detailed research and in 
this work we restrict ourselves to PageRank and CheiRank 
vectors. 
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Fig. 9. Spectrum of eigenvalues A of the Google matrix G of 
Wikipedia at different years. Red dots are core space eigen- 
values, blue dots are subspace eigenvalues and the full green 
curve shows the unit circle. The core space eigenvalues were 
calculated by the projected Arnoldi method with Arnoldi di- 
mensions ua ~ 6000. 



Finally we note that the fraction of isolated subspaces 
is very small for G matrix. It is increased approximately 
by a factor of order 10 for G* but still it remains very small 
compared to the networks of UK universities analyzed in 
[15] . This fact reflects a strong connectivity of network of 
Wikipedia articles. 



6 Discussion 

In this work we analyzed the time evolution of ranking of 
network of English Wikipedia articles. Our study demon- 
strates the stability of such statistical properties as PageR- 
ank and CheiRank probabilities, the article density distri- 
bution in PageRank-CheiRank plane during the period 
2007 - 2011. The analysis of human activities in different 
categories shows that PageRank gives main accent to pol- 
itics while the combined 2DRank gives more importance 
to arts. We find that with time the number of politicians 
in the top positions increases. Our analysis of ranking of 
universities shows that on average the global ranking of 
top universities goes to higher and higher positions. This 
clearly marks the growing importance of universities for 
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Fig. 10. Same as in Fig. [5]but for the spectrum of matrix G* 



the whole range of human activities and knowledge. We 
find that Wikipedia PageRank recovers 70 - 80 % of top 
10 universities from Shanghai ranking |23j . This confirms 
the reliability of Wikipedia ranking. 

We also find that the spectral structure of the Wikipedia 
Google matrix remains stable during the time period 2007 
-2011 and show that its arrow star structure reflects cer- 
tain features of small size unistochastic matrices. 
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7 Appendix 

The tables with all network parameters used in this work 
are given in the text of the paper. The notations used in 
the tables are: N is network size, N^ is the number of 
links, UA is the Arnoldi dimension used for the Arnoldi 
method for the core space eigenvalues, Nd is the number 
of invariant subspaces, dmax gives a maximal subspace di- 
mension, A'circ. notes number of eigenvalues on the unit 
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N 


Ne 


riA 


2003 


455436 


2033173 


6000 


2005 


1635882 


11569195 


6000 


2007 


2902764 


34776800 


6000 


2009 


3484341 


52846242 


6000 


200908 


3282257 


71012307 


6000 


2011 


3721339 


66454329 


6000 



Table 1. Parameters of all Wikipedia networks at different 
years considered in the paper. 





Ns 


iVd 


timax 


iVcirc. 


Nj 


2003 


15 


7 


3 


11 


7 


2003* 


940 


162 


60 


265 


163 


2005 


152 


97 


4 


121 


97 


2005* 


5966 


1455 


1997 


2205 


1458 


2007 


261 


150 


6 


209 


150 


2007* 


10234 


3557 


605 


5858 


3569 


2009 


285 


121 


8 


205 


121 


2009* 


11423 


4205 


134 


7646 


4221 


200908 


515 


255 


11 


381 


255 


200908* 


21198 


5355 


717 


8968 


5365 


2011 


323 


131 


8 


222 


131 


2011* 


14500 


4637 


1323 


8591 


4673 



Table 2. G and G* eigespectrum parameters for all Wikipedia 
networks, year marks spectrum of G, year with star marks 
spectrum of G*. 



circle with \Xi\ = 1, A^i notes number of unit eigenvalues 
with Aj = 1. We remark that Ns > TVcirc. > Ni > Nj. and 
Ns > dmax- The data for G are marked by the correspond- 
ing year of the time slot, the data for G* are marked by the 
year with a star. Links cleaning procedure eliminates all 
redirects (nodes with one outgoing link), this procedure is 
slightly different from the one used for the slot 200908 in 
[5]. All data sets and high resolution figures are available 
at the web page [26] . 
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